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I. INTRODUCTION 

The goal of using entangled electron pairs for the pro- 
cessing of quantum information poses a technological 
challenge that requires novel ideas on electron quantum 
transport. It has been proposed that a conventional 
superconductor is a natural source of entangled elec- 
trons which may be emitted into a normal metal through 
a properly designed interfaceii2i2ii£i£iL£i2ii2iii. At low 
temperatures and voltages, the electric current through 
a normal-superconducting (NS) interface is made exclu- 
sively of electron Cooper pairs whose internal singlet cor- 
relation may survive for some time in the context of 
the normal metal. The emission of two correlated elec- 
trons from a superconductor into a normal metal is of- 
ten described as the Andreev reflection-^ of an incident 
hole which is converted into an outgoing electron. The 
equivalence bet wee n the two pictures has been rigorously 
pro\ ed in Refs. l7l8lH^ . There the relation was established 
between the various quasiparticle scattering channels as 
these are referred to different choices of normal metal 
chemical potential, i.e. to different definitions of the vac- 
uum. When the reference chemical potential employed to 
label quasiparticle states in the normal metal is identical 
to the superconductor chemical potential (fiN — Ms), the 
number of Bogoliubov quasiparticles is conserved and the 
Andreev picture holds. If, on the contrary, ///v is chosen 
to be smaller than /j,g, quasiparticle number conservation 
is not guaranteed and spontaneous emission of two elec- 
trons through the SN interface becomes possible 8 . Trans- 
port calculations across an SN interface at low temper- 
ature and voltage which invoke an explicit two-electron 
picture have been presented in Refs. HHall4 



The need for spatial separation of the entangled beams 
has motivated the search for schemes that constrain (or 
at least allow) the two pair electrons to be emitted from 
different locations at the NS interface 1 -. In the conven- 
tional picture where quasiparticle scattering is unitary, 
that process is viewed as the absorption of a hole and 
its subsequent reemission as an electron from a distant 
point. Such a crossed (or nonlocal) Andreev reflection 
has been observed experimentalljii^iiSiil. 

The requirement of physical separation is a severe lim- 
itation in practice, since pairing correlations decay with 
distance. As a consequence, the current intensity of non- 
locally entangled electrons decreases with the distance 
r between the two emitting points. There is an ex- 
ponential decay on the scale of the superconductor co- 
herence length which reflects the short-range character 
of the superconductor pairing correlations 1 ^. A more 
important limitation in practice comes from the pref- 
actor, which, besides oscillating on the scale of the su- 
perconductor Fermi wavelength, decreases algebraically 
with distance. In the tunneling limit, and for a bal- 
listic 3D superconductor, the decay law is r~ 2 , if the 
tunneling matrix elements are assumed to be momen- 
tum independent 1 , or r -4 , if proper account is taken 
of the low-momentum hopping dependenc e) 8 ! 18 . Within 
the context of momentum-independent tunneling mod- 
els, the power law changes if the superconductor is low 
(d) dimensional^ 9 , or diffusiv o 5 i 20 , yielding r~ d+1 and 
r _1 , respectively. It remains to be investigated how that 
behavior changes when more realistic tunnel matrix ele- 
ments are employed 8 -* 1 " 8 " and when geometries other than 
planar or straight boundaries are considered. 

In this article we propose an experimental setup that 
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would guarantee long term separation of correlated elec- 
tron pairs without the shortcomings caused by the need 
to emit the pair electrons from distant points. The idea is 
to transmit both electrons through the same spatial region 
but inducing them to leave in different directions. In a 
ballistic normal metal such as a two-dimensional electron 
gas (2DEG), that divergent propagation guarantees the 
long term separation of the entangled electrons at dis- 
tances from the source much greater than the size of the 
source. 

To force the pair electrons to leave in different direc- 
tions, we propose to exploit the formation of resonances 
in a properly designed normal-superconductor interface. 
These could be one-electron (normal) resonances, such as 
those found in double-barrier structures^ (SININ struc- 
ture), or two-electron (Andreev) resonances such as the 
de Gennes - Saint-James resonances appearing in struc- 
tures with one barrier located on the normal metal side at 
some distance from the transmissive SN interface (SNIN 
structure )22iSi2ii2^. Those quasi-bound states have it in 
common that, in a perfect interface, they select the per- 
pendicular energy of the exiting electrons while ensuring 
the conservation of the momentum parallel to the inter- 
face. At low voltages and temperatures, this also deter- 
mines the parallel energy, given that the total energy of 
the current contributing electrons is constrained to lie 
close to the normal Fermi level. Altogether, this mecha- 
nism fixes the magnitude of the exit angle, since the par- 
allel momenta of the pair electrons are opposite to each 
other and both remain unchanged during transmission 
through the perfect interface. Thus the electron veloci- 
ties form a V-shaped beam centered around the perpen- 
dicular axis. 

The type of structures which are needed seems to 
be within the reach of current experimental expertise. 
In the last fifteen years, several groups have built a 
variety of hybrid superconductor-semiconductor (SSm) 
structureo 16 i 24 i 25 i 26 i 27 i 28 i 29 i 30 i 31 . More recently, some 
experimental groups 3 ^ 3 ^ have investigated transport 
through SSm structures where Sm is a 2DEG on a plane 
essentially perpendicular to the superconductor bound- 
ary. In such setups, the SN interface lies at the one- 
dimensional (ID) border of the two-dimensional (2D) 
ballistic metal. If two parallel straight-line barriers were 
drawn in that structure, one along the SN interface and 
another one at some distance within N, then the exper- 
imental scenario considered in this article would be re- 
produced. A three-dimensional (3D) version of the same 
structure, in which Sm would be 3D and the interface 
would be 2D, of the type reported in Ref. I2H would 
also produce divergent electron beams. These, however, 
would be emitted into a 3D semiconductor, where it may 
be more difficult to pattern suitable detectors. 

Once the two electrons propagate in the ballistic 
2DEG, their motion can be controlled by means of exist- 
ing techniques. For instance, they can be made to pass 
through properly located narrow apertures, such as those 
used in electron focusing experiments 3 -. For quantum 



information processing, their spin component in an arbi- 
trary direction could eventually be measured by using the 
Rashba effect* 3 ^ 7 , to rotate the spin before electrons en- 
ter the spin filter—. Then one could attempt to measure 
Bell inequalitiea 2 i 4 i 7 i 39 i 40 i 41 i 42 i 43 . Alternatively, one may 
measure electric current cross-correlations 9 * 2 £i 4 ^i^ii& to 
indirectly detect the presence of singlet spin correlations. 

In Section II, we describe the model we have adopted 
for our calculations. Two important features are the off- 
set between the conduction band minima and the differ- 
ence in the effective masses of S and Sm. Both effects 
have been analyzed by Mortensen et alM< in the context 
of SIN structures, with N a 3D semiconductor. In Section 
III, we focus on the linear regime and calculate the zero 
bias conductance using the multimode formula derived by 
Beenakker— . There we investigate the angular distribu- 
tion of the outgoing electron current and observe how it 
is indeed peaked around two symmetric directions. Sec- 
tion IV is devoted to the nonlinear regime 4 ^, where the 
voltage bias may be comparable to the superconductor 
gap. We find divergent beams again, this time with new 
features caused by the difference between the electron 
and hole wavelengths. By plotting the differential con- 
ductance, we relate our work to the previous literature 
on SN transport and note the presence of a reflcctionlcss 
tunneling zero bias peakSSiSL^S, as well as the existence 
of de Gennes - Saint-James resonances. In Section V, we 
discuss how the need to have a broad perfect interface, as 
required for parallel momentum conservation, can be rec- 
onciled with the interface finite size which is needed for 
the eventual spatial separation of the emerging beams. 
We conclude in Section VI. 



II. THE MODEL 

We wish to investigate the role of resonances in the 
angular distribution of the normal current in suitably de- 
signed SSm interfaces. A prototypical structure is shown 
in Fig. la, where the 2DEG forms an angle with the pla- 
nar boundaryof a superconductor, similar to the setup 
built in Ref. [M 

In the present analytical and numerical work we con- 
sider a semi-infinite ballistic 2DEG (hereafter also re- 
ferred to as N) lying in the half-plane x > 0. We assume 
a perfect interface, so that the one-electron potential is 
independent of y. Specifically, V(x) is taken of the form 

V(x) = -V e(-x) + HiSix) + H 2 S(x - L) . (1) 

Here, Vq accounts for the large difference between the 
widths of the S and N conduction bands. If Ep = 
h 2 kp/2m and E' F = h 2 kp/2m' are N and S Fermi ener- 
gies, respectively, one typically has E' F ~ Vq ^> Ep ^> A, 
where A is the zero-temperature superconducting gap. 
We assume that the bulk parameters change abruptly at 
x = 0. The structure contains two delta barriers, lo- 
cated at the SN interface and at a distance L from it 
within the N side. Their reflecting power is measured 
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FIG. 1: (Color online) (a) Schematic representation of the 
SN structure considered in this article. S is a conventional 
superconductor; N is a two-dimensional electron gas. Energy 
selection in the x direction by the resonance structure limits 
current flow to two divergent beams of entangled electrons, 
(b) One-electron potential profile in the x direction. Quasi- 
bound electron states form due to multiple reflection by two 
barriers of strength Z\ and Z2. 



by the dimensionless parameters Z\ and Z%, denned as 

1/2 

Z\ = Hi/h(vpv' F ) and Z2 = H2/hvp. The effective 
mass m, the Fermi wavevector /cf , and the Fermi velocity 
Vf are those of the normal 2DEG, while to', k' F , and v' F 
correspond to a conventional superconductor. 

It was shown in Ref . that the picture of two-electron 
emission and hole Andreev reflection are equivalent. For 
computational purposes, we employ here the standard 
Andreev picture whereby all quasiparticles have positive 
energy (e > 0), with the quasiparticle energy origin given 
by fj>s- However, in our discussion we will occasionally 
switch between the two images. An important feature 
is that the absence of a hole at e > in the Andreev 
scenario corresponds to the presence of an electron at 
— £ < in the two-electron picture^. 

In a transport context, the superconductor and normal 
metal chemical potentials differ by us — = eV, where 
V is the applied bias voltage. In the Andreev picture, one 
artificially takes [is as the reference chemical potential for 
labeling quasiparticles and the imbalance eV is accounted 
for by introducing an extra population of incoming holes 
with energies between and eV^^. 

An apparent shortcoming of the Andreev picture is 
that it does not show explicitly that the emitted elec- 
tron pairs are internally entangled. In this respect, we 
may note the following remarks: (i) the two-electron hop- 
ping matrix element vanishes when the spin state in the 
N side is a triplet 1 ; (ii) an analytical study of transport 
through a broad SN interface based on a two-electron 
tunneling picture 8 (with the final state explicitly entan- 



gled) gives results identical to those obtained within an 
Andreev description^; (iii) entanglement in the outgo- 
ing electron pairs has been explicitly proven in the gen- 
eral tunneling casei^; and (iv) transport across the SN 
structure is spin independent and thus must preserve the 
internal spin correlations of the emitted electron pair—. 
Moreover, using full counting statistics SamuelssonS has 
shown that current through an SN double-barrier struc- 
ture is carried by correlated electron pairs. 

To compute the current, we must sum over momenta 
parallel to the interface, which on the N side take values 
— kp < k y < kp. For the purposes of solving the one- 
electron scattering problem, we assume that the super- 
conductor is also two-dimensional. Due to the mismatch 
in effective masses, the perpendicular energy is not con- 
served (refraction) . The conserved quantum numbers are 
the parallel momentum [k y — k' y ) and the total energy 
(E x + E y = E' x + E' y -V = E, with E x £ E' x ). 

For a given k y , the energy available for perpendicular 
motion is E x — E — h 2 k y /2m, where E is the electron 
total energy. As a consequence, for each k y the picture 
depicted in Fig. lb holds provided that the (In is re- 
placed by an effective valued 



VN(k y ) = 1-J.n ~ h 2 k y /2m 



(2) 



which is matched to ns(k y ) = fis ~ ti 2 ky/2m\ with 
^s{k y ) — fiff{k y ) generally not equal to eV. 

Beenakker— has computed the SN zero bias con- 
ductance for an interface with many transverse modes. 
Mortensen et ai4£ have adapted the work of Ref. OH] 
to account for the full 3D motion through a perfect, 2D 
SSm interface, where the effective masses and the Fermi 
wavelengths of N and S may differ widely. Lesovik et 
alm^ have generalized the work of Refs. I48ll5ll to the 
nonlinear case where eV may be comparable to A. They 
have applied their results to structures displaying quasi- 
particle resonances. Here we combine the work of these 
previous three references. Specifically, we investigate the 
transport properties of an SN interface for arbitrary bias 
V between and A. We consider structures display- 
ing resonances due to multiple quasiparticle reflection, 
and allow for a large disparity between the S and N bulk 
properties. Most importantly, we calculate the angular 
distribution of the pair electron current emitted into the 
semiconductor. Another novel feature is that the semi- 
conductor we consider is a 2DEG whose plane forms an 
angle with the superconductor planar boundary, so that 
the SN interface is formed by a straight line. 



III. ZERO BIAS CONDUCTANCE 



The zero bias conductance is defined as 



G(0) ee lim dl/dV , 



(3) 
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FIG. 2: (Color online) Normalized angular distribution of the 
zero bias conductance for several values of the interbarrier 
distance L. The barrier strengths are (Z\, Z2) = (4, 2). Upper 
inset: schematic representation of the perpendicular potential 
profile. Lower inset: total conductance, normalized to half the 
maximum possible conductance, as a function of L. See main 
text for the meaning of other parameters. 



where / is the total current at voltage bias V. For an SN 
interface^, 



G(0) 



4e^ 



N 



Tr, 



(4) 



where {T„} are the eigenvalues of the one-electron trans- 
mission matrix through the normal state structure at to- 
tal energy E — /ijy — ^5 = /i, and N is the number of 
transverse channels available for propagation in the nor- 
mal electrode at energy fj,. For a perfect interface, the 
index v runs over the possible values of k y . Thus, when 
needed, we make the replacement — * (w/2ir) J dk y , 
where w — > 00 is the interface length. The minimum en- 
ergy required for propagation in mode v, referred to the 
bottom of the conduction band, is e v = h 2 k 2 /2m. 

In the linear regime, the total energy is restricted to 
be at fi. Therefore, the running value of k y determines 
the exit angle 



arctan (k y /k x ) 



since k x and k y must satisfy 



k 2 + ky 



kp 



Therefore, Eq. Q may be written as 

r 7r/2 



G(0) 



tt/2 



d9 G(0,0) 



(5) 
(6) 

(7) 



with G(0, 9) properly defined as the angular distribution 
of the zero bias conductance. 

In Fig. 2, we show G(0, 9) for several values of the in- 
terbarrier distance L, on a structure with potential bar- 
riers of strength Z\ = 4 and Z2 = 2 located at x = 



and L, respectively. It is divided by Ae 2 w/h\p (Xp be- 
ing the N Fermi wavelength) , which is half the maximum 
possible value of G(0) (obtained when T v = 1 for all v). 

The semiconductor conduction band width is taken 
Ep =ijx 100 K. The ratios between the Fermi wavevec- 
tors and Fermi velocities in N and S are, respectively, 
rfc = kp/k'p = 0.007 and r v = vp/v' F = 0.1 (GaAs val- 
ues). The presence of quasi-bound states located between 
the two barriers yields a structure of resonance peaks in 
the one-electron transmission probability T v as a function 
of v. We also note that the small value of r^ will cause 
important internal reflection of the electrons within the 
superconductor. As a result, only S electrons very close 
to normal incidence will have a chance to be transmitted 
into N. Once in N, they may leave with much larger an- 
gles. Specifically, if 9' is the angle on the S side, one has 
sin#' = r^sin^ (Snell law). For the parameters consid- 
ered in this article, only electrons arriving from S within 
A6*'/2 = arcsin(r-fc) ~ 0.4 degrees of normal incidence are 
transmitted through the normal-state structure. 

As L increases, the position of the resonant levels is 
lowered. In Fig. 2, the values of L are chosen such that 
only the lowest resonant level plays a role. This allows us 
to investigate the effect of a resonant level at perpendic- 
ular energy (on the N side) E x = En < M> which appears 
as a peak in T„ as a function of v. This occurs for v = Vr 
satisfying 



H - 



= E h 



(8) 



For the shortest interbarrier distance displayed (L = 
23 nm), the structure of G(0, 9) begins to reveal the pres- 
ence of a resonance just below Ep. The trend towards a 
bifurcation of the conductance angular distribution be- 
comes clearer for larger values of L. As discussed before, 
the presence of a sharp resonance only permits the trans- 
mission of electrons with perpendicular energy E x close 
to En. This fixes the value of k x at k x = kf> and, with 
it, the magnitude of the exit angle 



9 R = arctan ^y/(k F /k R ) 2 



1 



(9) 



For a given linewidth T of the one-electron resonance, the 
corresponding spread of the angular distribution is 



A9 



r 



E F sm{29 R ) ' 



(10) 



Thus, the angular width has a minimum at 9 R = 7r/4, as 
in fact revealed by the narrower spikes in Fig. 2. 

The lower-right inset of Fig. 2 shows the total con- 
ductance [see Eq. 0] as a function of the interbarrier 
distance. It is normalized to half its maximum possible 
value. For small L, the lowest resonance lies at E R > /i, 
which blocks current flow. As L is increased, E R de- 
creases and the lowest resonance becomes available for 
transport (E R < //). Then G(0) shows a rapid increase 
followed by a decaying tail. The effect is so marked that, 
if we attempt to plot G(0, 9) for e.g. L = 22 nm (just 
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FIG. 3: (Color online) Same as Fig. H for (Z u Z 2 ) = (2, 2). 



SNIN 



L 

G(0)/[(4e 2 /h)(w/A F )] 


0.3 








" J 










\\\ \ 

\ N v V 



20 40 60 

L{nm) 




////// T-=1 

'\t//y r ^ - 00 " 



00K 



007 
r=0.1 

Z,=2 



-0.60 0.00 0.60 

FIG. 4: (Color online) Same as Fig. H for (Zi, Z 2 ) = (0, 2). 



below the smallest shown value), the resulting curve is 
invisible on the scale of Fig. 2. As L increases further, 
a second resonance becomes available for transmission 
and the wide spikes due to the the first resonance coexist 
with the new, more centered lobes which in turn tend to 
bifurcate as L increases even more (not shown). 

The decay of G(0) for L > Lr (where Lr is the in- 
terbarrier distance at which Er = /i) goes like L^ 1 / 2 
because it reflects the ID nature of the transverse den- 
sity of states. This can be proved by noting that Eq. Q 
can be written as 



4e 2 x - 



(ii) 



where A v = T 2 /(2 — T u ) 2 is the probability for Andreev 
reflection in mode v at total energy /i, which corresponds 
to quasiparticle energy e = 0. Because of the normal 
resonance, both T v and A v are strongly peaked around 
the value of vr satisfying JSJ. Thus we may approximate 
A v ~ aS(fi — e v — Er), where a is an appropriate weight. 
Then G(0) becomes 



G(0) ~ ^-aD(fi - E R ) 
h 



(12) 



where D(e) = ^„ S(e — e v ) is the transverse density of 
states. On this energy scale, Er is a smooth function of 
L, so that it can be approximated as Er ~ /j, — b(L — LR), 
with b > 0. Then Eq. (Tg) yields G(0) oc D(b(L-L R )) ~ 
(L — Lfl) -1 / 2 , as observed in the inset of Fig. 2. Such 
a manifestation of the transverse density of states in the 
total transport properties is characteristic of structures 
which select the energy in the propagation perpendicu- 
lar to the plane of the heterostructure^. The foregoing 
argument allows us to predict that, for a 3D structure, 
the total conductance will display steps as a function of 
L, since then D(e) will be constant (not shown). 

Figs. 3 and 4 show G(0, 9) for setups identical to that 
of Fig. 2, except for Z\ taking values 2 and 0, respec- 
tively, Zi remaining fixed at 2. The building of SSm 



interfaces with small Z\ seems feasible with the doping 
techniques implemented in Refs. I25ll2^l30l 

As in Fig. 2, the electron flow is channeled through 
well-defined resonances in the x direction, again giving 
rise to divergent beams in the N electrode. At first sight it 
may seem surprising that for Z\ = one still finds peaks 
in the angular distribution, since they reveal a structure 
in the transmission T v that is not expected from a single 
barrier of strength Z^. However, when Z\ = 0, there is 
still some normal reflection at x = due to the large 
mismatches E' F 3> Ep and m! 3> m. In fact, on quite 
general grounds, one has T v — > as e u — > Ep (equiva- 
lent to k x — > 0), even if Z\ =0. This trend is revealed 
by the decreasing length of the spikes for increasing 
(decreasing k x ). 

In Fig. 3, fi stays slightly above Er for L = 23 nm. 
The details of reflection at the interface cause some shift 
in the detailed position of the resonances. For Z\ = 
(Fig. 4), the resonant level Er at that particular in- 
terbarrier distance is exactly at fi, as revealed by the 
absence of splitting in G(0,8). If, by decreasing L, Er 
were taken considerably above //, then the forward lobe 
of Fig. 4 would be sharply reduced. This general prop- 
erty was already noted in the discussion of Fig. 2 and its 
inset. 



IV. NONLINEAR TRANSPORT: SPECTRAL 
CONDUCTANCE 

We have seen that, in the zero bias limit, the peaks in 
the angular distribution directly reflect the structure of 
(normal) resonances in T u as a function of v, since this 
determines G(0) through Eq. As V becomes nonzero 
and comparable to A, new resonances appear which are 
a direct manifestation of Andreev reflection occurring at 
nonzero quasiparticle energies. Such Andreev resonances 
have been discussed, for instance, in Refs. I23u25ll49l Be- 
low we present a brief description that suits our present 
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needs and which complements the discussion given by 
Lesovik et al~. 

We restrict our study to the case < \eV\ < A. As 
in Ref. 0, we focus for simplicity on the spectral con- 



G(e,V) 
9u(e,V) 



Here, g v (e, V) is the Andreev reflection probability for 
a quasiparticle of energy e incoming in mode v, with 
e < \eV\. It is determined by T v (e), which is de- 
fined as the transmission probability for an electron inci- 
dent from the N side on the normal structure (i.e. with 
A = 0) in transverse mode v with total energy \i$ + e, 
R v {s) = 1 — T v {e), i?(e) = arccos(e/A), and yv(e) is 
the phase of the reflection amplitude for an electron im- 
pinging from the S side on the normal structure. The 
latter depends on e through the phases acquired upon 
reflection on each barrier (usually negligible) and, more 
importantly, through the optical path between the two 
barriers k u (e)L, where 

fc„(±e) = [2m(E F + eV - e v ± e)] 1/2 /h 

~ mvp v /h+ (eV ±£)/hvp u . (15) 

1/2 

Here, vp v — [2 [Ep — e„) /m] ' is the perpendicular ve- 
locity for a Fermi electron in mode v on the N side (note 
that Ep+eV is the energy difference between the S chem- 
ical potential and the bottom of the N conduction band). 
We notice the symmetry g(e, V) = g{—£, V) and the fact 
that, through l|15|l. the transmission T v (e) does depend 
on V . In practice, we are only interested in the case 
e = eV. Thus, hereafter we refer to both G and g as 
functions of a single argument e which is to be identified 
with eV in the sense indicated in Eqs. I|13(l and 1)14(1. 

The structure of the angular distribution of the con- 
ductance reflects that of g v as a function of v, which 
generally reveals a complex and rich behavior, since 
it is determined by the combined role of the product 
T v (e)T v (—e) and the cosine term in l|14|). Below we dis- 
cuss some general trends. 

First we note that g v (0) = T 2 (0)/[2 - T„(0)} 2 , with 
T y (0) computed for eV = 0, which is consistent with Eq. 
Q. If the one-electron (normal) resonance occurs at a 
perpendicular energy E x — Er satisfying hn — Ep < 
Er < fJ,s, for |e| < A there is always a transverse mode 
v(e) for which 

Ms - e»(e) +e = E R , (16) 

i.e. such that T v {e) presents a peak at v = v(e) as a func- 
tion of v, with maximum value Tq (normal resonance). 



ductance G(e, V), i.e. we neglect the contribution to the 
total differential conductance coming from the derivative 
with respect to V of G(e, V) itself. From Ref. we 
note that, for < |e| < A, 



(13) 
(14) 



In a symmetric structure, Tq = 1. 

As a function of v, the phases ip u (ie) undergo an 
abrupt change near i/(±e), so that the cosine term goes 
quickly through two maxima, in v(e) and v{— e), none 
of which necessarily reaches unity. These maxima co- 
incide in general with the peaks of T v (e) and T u (—e). 
From (|14|) . this translates into pairs of close lying peaks 
in the conductance angular distribution. We have ob- 
served that the above tendency is typically present for 
all intermediate values of e (as compared with A) for 
{Z\,Z2) = (4,2) and (0,2). Now we describe another 
aspect of the peak formation mechanism that is relevant 
for e not much smaller than A in the structure (0,2). 
We note that it is compatible with the trend discussed 
above. 

Andreev resonances are characteristically given by the 
condition^ 

cos [¥>„(£)-?>„(-£)- 20(e)] = 1 . (17) 

If we recall that e is to be identified eventually with eV, 
and that through (|15fl f 1/ {e) does also depend on V, we 
may state that, for a continuous range of voltages V, 
there is always at least a value of v = v{V) satisfying 
(|17|l . As defined in (|16fl . v({)) is also a function of volt- 
age, since /is = fxpf + eV with /x^r fixed. Alternatively, 
one may take us as fixed and /zjv dependent on voltage; 
then, i/(0) is independent of V. In both scenarios (and, 
conceivably, in intermediate ones), there is a discrete set 
of values {V n } for which the two transverses modes coin- 
cide, i.e. , for which D(V n ) — v(0). 

We note on the other hand that, for e = 0, <|16[) may 
also be regarded as the maximum condition for T y ( )(e) 
viewed as a function of e with its maximum lying at 
e = (i.e. at total energy us)- Thus we can assert that 
^i>(o)( £ ) — T u (q)(— s) within a range of e values, which 
may include e n = eV n . Noting that the Andreev reso- 
nance condition 117|) is symmetric in e, we conclude from 
iJUJl that 

01/(0) fen) = 1 . (18) 

even if T l/ ( )(£ n ) is not unity. This maximum value of 
the conductance per mode (which is 2 in units of 2e 2 /h; 



T v {e)T v {-e) 



1 + R u {e)R v {-e) - 2 [R w {e)R v {-e)] 1/2 cos [^(e) - p„(-e) - 20(e)] 

I 
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FIG. 5: (Color online) Normalized angular distribution of the 
spectral conductance evaluated at e = 0.4 meV, for several 
values of the interbarrier distance L. The superconductor gap 
is A = 1 meV. Barrier strengths are (Z\,Z'i) = (4,2). Left 
inset: schematic representation of the perpendicular potential 
profile. Right inset: total current, integrated over angles and 
energies (up to A), normalized to half its maximum possible 
value, as a function of L. 



see Ref. |57J for a discussion) is consistent with the results 
reported in the single-mode study of Ref. |23|. Therefore, 
at voltages {V n } the total transmission (summed over 
v) receives a strong contribution from v = ^(0) and its 
vicinity. This behavior tends to generate peaks in the 
total spectral conductance G(e) at or near the values 
£n = eV n defined above. 

The conclusion is that the sharpest resonances nucleate 
at angles near normal resonances {v(e) is typically close 
to vr, since |e| < A <C Ep). This happens for all energies 
e. However, as explained above, some energies e benefit 
more efficiently from the resonance (in the sense that 
g v (s) displays higher maximum values as a function of 
v) and thus give rise to peaks in G(e) when integrated 
over angles. 

Now we may argue like in Section III. Whenever [is > 
En, there is a low-lying transverse mode v satisfying (11611 . 
Then we expect to have a strong peak in the angular 
distribution of the spectral conductance, G(e,0), which is 
defined to yield 



r ir/2 

G(e) = I d6 G(e, ( 

-7T/2 



(19) 



Figs. 5-7 show the normalized value of G(e, 9) for struc- 
tures with (Zi,^) = (4,2) and (0,2), the former being 
considered for two different combination of e and A. As 
L increases, the value of Er decreases and sinks below 
fis- This generates maxima in the angular distribution 
in the manner discussed above. 

At zero temperature, and for eV > 0, G(e) can be un- 
derstood as the contribution to the total current stem- 
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FIG. 6: (Color online) Same as Fig. [S] for A = 0.1 meV and 
e = 0.099 meV. 
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FIG. 7: (Color online) Same as Fig. for (Z U Z 2 ) = (0,2) 
and e = 0.45 meV. 



ming from electron pairs emitted into the normal metal 
with total energies ± £■ The two electrons leaving the 
superconductor have identical \k x \ and slightly different 
total energy (see below). Thus they do not point exactly 
in the same direction, i.e. the V which they form upon 
emission is not exactly centered around the normal axis. 
By symmetry, for each pair in which e.g. the upper elec- 
tron is emitted towards the right (and the lower one to 
the left), there is another pair solution in which the up- 
per electron travels to the left (and the lower one to the 
right). When plotting the total differential conductance, 
the two asymmetric Vs appear as a single V whose lobes 
are double peaked. 

We note here that, in the contribution to G(e) as de- 
fined in Eqs. l(T3|l and ifHJl. T v {— e) is identical to the T v 
appearing in the zero voltage limit discussed in the pre- 
vious section [see Eq. lfT5|) ]. i.e. v{— e) = as defined 
in JSJl, if we identify (In = M- This implies that, in the 
double-peaked lobes, the inner peak points in the same 
direction as the single-peaked lobe of the linear (V = 0) 
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FIG. 8: (Color online) Subgap spectral conductance (inte- 
grated over angles), normalized to half its maximum possible 
value, for four values of the interbarrier distance L, two values 
of the superconductor gap A, and two values of the strength 
Z\ of the barrier located right at the SN interface. 



limit, a result which is independent of the sign of eV. 
The fact that the coincidence occurs at the inner peak 
can be understood by noting that, since e = eV , we have 
AV(e)(e) = fe v (_ e )(— e), while e u{e) = e u{ _ e) + 2eV. Thus, 
at a given s, peaks in the angular distribution occur at 
v{e) and v{—e). Both have the same perpendicular mo- 
mentum, but the latter has lower parallel kinetic energy. 

The fact observed in Figs. 5-7 that the inner peak 
displays a larger current density is due to the asymmetric 
character of the peaks in T„(±£) as a function of v (or the 
angle 9), which ultimately reflects the greater efficiency 
with which close-to-normal emission electrons contribute 
to the electric current. 

The insets of Figs. 5-7 show the total current (inte- 
grated over 9 and e) as a function of L. As for the zero 
bias conductance, they reveal a succession of maxima fol- 
lowed by an inverse square root decay law that mirrors 
the transverse density of states (see discussion in the pre- 
vious section). 

Fig. 8 shows the total spectral conductance for volt- 
ages below the gap. This type of curves has been the 
object of preferential attention in the previous literature 
on NS transport. By presenting them here, we make 
connection with that preexisting body of knowledge, in 



particular with the experimental and theoretical works 
of Refs. and lifi respectively. The forthcoming re- 
marks are intended to complement that discussion and 
to provide a self-contained, unified picture of the work 
presented here. 

The asymmetry in G{eV) is due to the finite normal 
bandwidth. For the results plotted in Fig. 8, the voltage 
V varies as /is varies with /i^r fixed. From Fig. 1, it 
is clear that raising /is is not equivalent to lowering it. 
Asymmetric curves are displayed in Ref. 03 and have 
been discussed in Ref. (see also references therein). 
In what follows, we focus on the behavior for eV > 0. 

Both in Figs. 8a and 8b we present two groups of 
curves, corresponding to a small and a large gap. The 
barrier parameters of Fig. 8a are the same as those of 
Figs. 5 and 6 , namely (Zx,Z 2 ) = (4,2). Although 
Figs. 5 and 6 already exhibit Andreev features such as 
the double-peaked lobes in G(e, 0), these are washed out 
when the angular variable is integrated to yield the total 
spectral conductance G(e = eV), as shown by the single 
peaked curves obtained for the same value of the gap as 
in Fig. 5 (A — 1 meV), or by the absence of peaks for 
the parameters of Fig. 6 (A = 0.1 meV). The curves for 
A = 1 meV display a clear zero bias conductance peak 
(ZBCP) whose height is determined by the structure nor- 
mal properties [see Eq. Q]. As e increases above zero, 
both electron and holes (or both the upper and lower 
energy emitted electrons) may benefit from the low-lying 
normal resonance (En < /is) as long as e < T, where T is 
the linewidth of the normal resonance. When e > T, it is 
not possible to channel both electrons through the same 
resonance and the contribution to the conductances de- 
creases. On closer inspection, one finds that the width of 
the ZBCP is indeed determined by the normal resonance 
width, but not by that appearing in the perpendicular 
transmission T v (e) (viewed as a function of e). Rather, 
it essentially mirrors the width of the numerator in Eq. 
(fHj) . This is the product Ty(e)T v (— e) evaluated at vr 
and viewed also as a function of e, i.e. for electrons leav- 
ing in the direction of maximum current flow (at exit 
angle 9 = 9r). This is reminiscent of the result stating 
that, when Zi is replaced by a disordered normal metal, 
the width of the ZBCP is of the order of the Thouless 
energy^. 

A general property of SN interfaces with a single bar- 
rier right at the interface is that Andreev reflection prob- 
ability tends to unity as |e| — > A£i. However, we find that 
this is generally not the case for a double barrier inter- 
face. For Z\ = 0, we do notice that sharp peaks in G(e) 
form just below the gap for some values of L, so close 
to it that they can be observed only through a magnifi- 
cation of Fig. 8. Due to this tendency to acquire large 
values near the gap, G(e) goes through a minimum at 
finite e if the width of the ZBCP is smaller than the gap. 
This is the case shown in Fig. 8a for A = 1 meV. For a 
smaller gap (A = 0.1 meV), the value of G(0) remains 
unchanged but there is no room for G(e) to display a 
minimum between and A. 
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Being Z\ — more transmissive (although not entirely, 
because of the reflection at the potential step; see Section 
III), Fig. 8b displays Andreev resonance features that do 
survive upon integration over angles. For A = 1 meV 
and L = 23 nm, one observes a peak at finite energies 
that adds to the overall ZBCP. As L increases, the inner 
Andreev peak evolves towards zero energy. At larger dis- 
tances (L = 36 nm), the lowest Andreev resonance can 
only be hinted at as a shoulder in the plot for A = 0.1 
meV. We also note that, for L = 24 and 26 nm, a sec- 
ond Andreev resonance becomes visible close to the gap 
edge. However, due to the involved interplay between 
the transmission probabilities and the cosine term ap- 
pearing in Eq. JOJl, this second peak does not appear to 
follow a simple monotonic trend. In fact, for A = 1 meV, 
the second resonance is no longer observable because it 
evolves towards a sharp peak just below the gap. 



V. DISCUSSION 

So far we have assumed that the SN interface is in- 
finitely long (w — > oo). This has allowed us to treat k y as 
a continuous, conserved quantum number, which consid- 
erably simplifies the transport calculation. Of course, the 
idea of an infinite interface is at odds with the primary 
motivation of our work, which is to propose a method 
to spatially separate mutually entangled electron beams. 
Below we argue that, fortunately, only a moderately long 
interface is needed in practice. 

For simplicity, we focus our discussion on the low volt- 
age limit, where the total energy can be assumed to be 
sharply defined. Then the width A9 of the angular dis- 
tribution is due only to the uncertainty in the parallel 
momentum Ak y . This in turn is closely connected to 
Ak x through the relation k x Ak x = k y Ak y , since total 
energy uncertainty is zero. There are two contributions 
to the momentum uncertainty: the nonzero width of the 
resonance in the perpendicular transmission and the fi- 
nite length of the SN interface. Thus we may estimate 



Afc„ 



mT 

fi? kit 



1 

w 



This translates into an angular width 



AO 



1 



Ef$oi(26r) kpwcosOR 



(20) 



(21) 



The actual angular width of G(0, 6) is actually a bit 
smaller, since the present estimate is based on one- 
electron considerations, while the relevant angular dis- 
tribution is determined by Eq. Q}. We neglect this dif- 
ference for the present simple estimates. 

Eq. 1)21(1 contains two contributions. The first term is 
determined by the normal resonance and is responsible 
for the width of the angular distributions plotted in Figs. 
2-4 (with w — > oo). Our main concern here is that the 
second contribution, that which stems from the finiteness 
of the aperture, does not contribute significantly. 



A strict criterion may be that the interface finite 
length should not modify the intrinsic angular width 
(Hvf sin 9r/w <C r), which everywhere has been assumed 
to be small enough to allow for narrow divergent beams. 
A more lenient criterion is that, regardless of the spe- 
cific value of r, the finite aperture should not generate 
an excessively broad angular distribution. For typical 
cases this amounts to requiring kpw 3> 1 (for a discus- 
sion see Fig. 5 in Ref. @). For the bandwidth which we 
have assumed (Ep/kB = 100 K) and an effective mass 
of m = (rk/r v )m' — 0.07m e , where m e is the bare elec- 
tron mass, we have Xf = 2it /kp ~ 50 nm. So apertures 
greater than a few hundred nanometers seem desirable to 
keep the angular uncertainty within acceptable bounds. 

Another source of angular spreading is interface rough- 
ness, with a characteristic length scale I. However, 
it should not pose a fundamental problem as long as 
I ^> Xf, so that a structure of intermediate width could 
be designed satisfying I 3> w 3> Xf- 

For the difference in velocity direction to translate into 
spatial separation, it is necessary that the spin detectors 
are placed sufficiently away from the electron-emitting 
SN interface. Of course, the needed distance depends also 
on the exit angle 9r. For a convenient value of Or ~ 7r/4, 
simple geometrical considerations suggest that, unsur- 
prisingly, the distance d from the detector to the center of 
the SN interface must be greater than its width w. Since 
elastic mean free paths in a 2DEG can be made as high as 
l e ^100 /im, there seems to be potentially ample room for 
building structures satisfying Xf -C w <C d <C l e - Such 
devices would display well-defined divergent current lobes 
which could be detected (and, eventually, manipulated) 
at separate locations before the directional focusing is 
significantly reduced by elastic scattering. 



VI. CONCLUSIONS 

We have investigated theoretically the possibility of 
creating hybrid normal-superconductor structures where 
the two electrons previously forming a Cooper pair in the 
superconductor are sent into different directions within 
the normal metal. The central idea relies on the design 
of a structure that is transparent only to electrons with 
perpendicular energy within a narrow range of a reso- 
nant level. Since the total energy lies close to the Fermi 
level, such a filtering of the electron perpendicular energy 
translates into exit angle selection. 

Electrons from a conventional superconductor are 
known to be correlated in such a way that electrons mov- 
ing at similar speeds in opposite directions tend to have 
opposite spin. At low temperatures and voltages, elec- 
tron flow from the superconductor to the normal metal is 
entirely due the transmission of correlated electron pairs. 
These have both opposite spin and opposite parallel (to 
the interface) momentum, while possessing the same to- 
tal energy. If the exit angle is selected by filtering the per- 
pendicular momentum, the current in the normal metal 
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is formed by two narrow, mutually singlet entangled elec- 
tron beams which point in different directions and which 
spatially separate from each other at distances from the 
source much greater than the width of the source. 

The trick of exit angle selection is intended to facil- 
itate a neat observation of nonlocal entanglement be- 
tween electron beams, and this article has been devoted 
to proposing a specific implementation of that idea. One 
cannot help noting, however, that such a selection of 
the outgoing direction might not be totally essential. If 
we content ourselves with measuring anticorrelated low- 
energy spin fluctuations over mesoscopic length scales, it 
may just be sufficient to place the two spin detectors sym- 
metrically around the interface at a sufficient distance 
and angle, very much like in the setup of Fig. la but 
with a conventional, non angle-selecting SN tunnel inter- 
face. If their motion between the emitter and the detec- 
tor is ballistic, electrons arriving at each detector have, 
on average, opposite parallel momentum and opposite 
spin (angular anticorrelation has been explicitly shown 
in Ref. [U for a broad perfect interface) . The boundaries 
of the 2DEG might conceivably be designed to optimize 
such correlations. The outcome is that electrons arriving 



at each detector will exhibit a degree of nonlocal spin- 
singlet correlations that could be measured. 

Altogether, we conclude that a ballistic two- 
dimensional electron gas provides an ideal scenario to 
probe nonlocal entanglement between electrons emitted 
from a distant, finite-size interface with a superconduc- 
tor. If that interface is formed by a hybrid structure that 
selects the perpendicular energy and thus the magnitude 
of the electron exit angle, nonlocal spin correlations will 
be clearly observed if the outgoing beams are directed 
towards suitably placed detectors. 
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